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Abstract. We study the Partial Degree Bounded Edge Packing (PDBEP) 
problem introduced in [5 > by Zhang. They have shown that this prob- 
lem is NP-Hard even for uniform degree constraint. They also presented 
approximation algorithms for the case when all the vertices have degree 
constraint of 1 and 2 with approximation ratio of 2 and 32/11 respec- 
tively. In this work we study general degree constraint case (arbitrary 
degree constraint for each vertex) and present two combinatorial approx- 
imation algorithms with approximation factors 4 and 2. We also study 
integer program based solution and present an iterative rounding algo- 
rithm with approximation factor 3/(1 — e) 2 for any positive e. Next we 
study the same problem with weighted edges. In this case we present an 
O(logn) approximation algorithm. Zhang [5] has given an exact 0(n 2 ) 
complexity algorithm for trees in case of uniform degree constraint. We 
improve their result by giving 0(n- log n) complexity exact algorithm for 
trees with general degree constraint. 
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1 Introduction 

The partial degree bounded edge packing problem (PDBEP) is described as fol- 
lows: Given a graph G = (V,E) and degree-bound function c : V — > N, com- 
pute a maximum cardinality set E' C E which satisfies the degree condition: 
(d' u < c u ) V (d' v < c v ) for each e = (u,v) E E'. Here d' x denotes the degree of 
vertex x in the graph G' — (V,E'). Without loss of generality, we will assume 
that c v < d v for all v S V where d v denotes the degree of v in G. 

In the weighted version of the problem edges are assigned non-negative 
weights and we want to compute a set of edges E' with maximum cumulative 
weight subject to the degree condition described above. 

In [5] , the PDBEP problem was motivated by an application in binary string 
representation. It was shown there that the maximum expressible independent 
subset (MEIS) problem on 2-regular set can be reduced to PDBEP problem 
with uniform constraint c = 2. The PDBEP problem finds another interesting 
application in resource allocation. Given n types of resources and m jobs, each 



job needs two types of resources. A job j, which requires resources u and v, can 
be accomplished if u is not shared by more than c u jobs or v is not shared by 
more than c v jobs. Interpreting the resources as the vertices of the input graph 
and the jobs as edges, the PDBEP problem is to compute the maximum number 
of jobs that can be accomplished. 

1.1 Related Work 

The decision problem of edge packing when there is a uniform degree constraint 
of 1 is a parametric dual of the Dominating Set (DS) problem. The parametric 
dual means that for graph G — (V, E), a k sized dominating set implies a \V\ — k 
sized edge packing, and vice versa. The parametric dual of DS was studied in 
[3J. Further, the dual was well studied under the framework of parameterized 
complexity by Dehne, Fellows, Fernau, Prieto and Rosamond in pQ. 

Recently Peng Zhang [5] showed that the PDBEP problem with uniform 
degree constraint (c v = k for all v) is NP-hard even for k = 1 for general graphs. 
They gave approximation algorithms for the PDBEP problem under uniform 
degree constraints of k = 1 and k = 2 with approximation factors 2 and 32/11 
respectively. They showed that PDBEP on trees with uniform degree constraint 
can be solved in 0{n 2 ) time. 

1.2 Our Contribution 

We propose three different approximation algorithms for the problem with gen- 
eral degree constraints (i.e., for arbitrary function c). Two of these algorithms 
are combinatorial in nature and their approximation ratios are 4 and 2. We then 
show that the relaxation of the natural integer program for this problem has a 
large integrality gap. Then we propose an "approximate" integer program which 
is a Lagrangian-like relaxation of the original IP, and show that any a approx- 
imation of this IP is a 2a/ (1 — e) approximation of the PDBEP problem for 
any e > 0. We then present a 1.5/(1 — e) approximation iterative rounding [5] 
algorithm for the new integer program. Although this only leads to a 3/(1 — e) 2 
factor approximation, we think that this Lagrangian like relaxation is an im- 
portant contribution and hope that this method can be applied to some other 
problems to get significantly better results especially in cases where the natural 
IP has a large integrality gap. 

The results detailed above are significantly improved results over the 2 and 
32/11 approximations in [5] which are applicable to constant function cases c = 1 
and c = 2 respectively. 

Next we consider the PDBEP problem with general degree constraint for 
edge-weighted graphs. In this case we present a combinatorial approximation 
algorithm with approximation factor of 2 + 2 logn. 

Finally we present an exact algorithm for unweighted trees with general de- 
gree constraint function. The time complexity of this algorithm is 0(n log n). 
This is an improvement over the 0(n 2 ) algorithm in [5] which is applicable to 
only a constant degree constraint function. 



2 Approximation Algorithms for the unweighted case 

The optimum solution of a PDBEP problem can be bounded as follows. 

Lemma 1. Let G = (V,E) be a graph with degree-bound function c : V — > N. 
Then the optimal solution of PDBEP can have at most Ylvev ° v e dges- 

Proof. Let E'cSbea solution of PDBEP. Let U = {v e V\d' v < c v }. Then 
from the degree condition we see that U is a vertex cover in the graph (V, E'). 
Hence \E'\ < Eue(7 c « ^ J2vev c ^- D 

2.1 Edge Addition based Algorithm 

Consider any maximal solution Y C E, i.e., 7U{e} is not a solution for e € E\Y. 
Let dy{x) denote the degree of a vertex x in the graph (V, Y). Partition the vertex 
set into sets: A — {v\dy{v) < c v }, B — {v\dy(v) = c v }, and C = {w|dy(w) > c v }. 
Observe that every edge of the set E\Y which is incident on a vertex in A, has 
its other vertex in B. Hence for any «i,a2 £ A the E\Y edges incident on a\ are 
all distinct from those incident on et 2 . Construct another edge set Z containing 
any c v — dy(v) edges incident on v for each v G A. Observe that Z also satisfies 
degree constraint. Output the larger of Y and Z. See Algorithm [TJ We have the 
following result about the correctness. 

Lemma 2. The Algorithm^ outputs a set which satisfies the degree constraint. 

Consider the set Y U Z. In this set the degree of each vertex is not less than 
its degree-bound. Hence the cardinality of the output of the algorithm is at least 
^2 V c„/4. From Lemma [T] the approximation ratio is bounded by 4. 

Theorem 1. Algorithm]]} has approximation factor 4. 

2.2 Edge Deletion based Algorithm 

The second algorithm, Algorithm^ for PDBEP is based on elimination of edges 
from the edge set. Starting with the input edge set E, iteratively we delete the 
edges in violation, i.e., in each iteration one edge (u,v) is deleted if the current 
degree of u is greater than c u and that of v is greater than c v . The surviving 
edge set Y is the result of the algorithm. Clearly Y satisfies the degree condition. 
Also observe that dy(v) > c v for all v £ V. Hence \Y\ > c u/2- From Lemma 
E \Y\ > OPT/2. 

Theorem 2. Algorithm^ has approximation ratio 2. 

2.3 LP based Algorithm 

In this section we explore a linear programming based approach to design an 
approximation algorithm for PDBEP. 



Data: A connected graph G = (V, E) and a function c : V — > N such that 

c v < d(v) for each vertex v. 
Result: Approximation for the largest subset of E which satisfies the 
degree-condition. 

Y :=0; 

for e € E do 

if Y U {e} satisfies the degree- condition then 
| Y:=YU{e}; 

end 
end 

Compute A := {v € V\dy{v) < c}; 
for v £ A do 

j Select arbitrary c v — dy (w) edges incident on v in _E \ Y and insert into Z; 
end 

if |y| > \Z\ then 

! return y"; 
else 

i return Z; 
end 

Algorithm 1: Edge Addition Based Algorithm 



Data: A connected graph G = (V, E) and a function c : V — > N such that c„ is 

the degree bound for vertex v. 
Result: Approximation for the largest subset of E which satisfies the 
degree-condition. 

Y:= E; 

for e = (u, v) € Y do 

if dy{u) > c u and dy(v) > c v then 
I Y^Y\{e}; 

end 
end 

return Y; 

Algorithm 2: Edge Deletion Based Algorithm 



The Integer Program Following is the natural IP formulation of the problem: 



i' = Ve 

e<EE 

He < Xu + %v Ve G E 

^2 He < c v x v +d v (l- x v ) Vv <EV 

eeS(v) 

x v G {0, 1} Vu G V 
y e G {0, 1} Ve G E 

The solution computed by the program is E' = {e\y e — 1}. The linear pro- 
gramming relaxation of the above integer program will be referred to as LP1. 

Lemma 3. The integrality gap of LP1 is fi{n) where n is the number of vertices 
in the graph. 

Proof. Consider the following instance of the problem. Let G be a complete graph 
on n vertices {vq, v%, . . . , v n -i} and the degree constraint be c v = 1 Vv G V . We 
now construct a feasible fractional solution of LP1 as follows. Let x v = 0.5 for all 
v and y e = 1 for all e = (i>i, Vj) where j is in the interval (i — \ n/A\ (mod n), i + 
[n/4j(mod n)). The value of the objective function for this solution is at least 
(n — l) 2 /4. On the other hand, from Lemma[l] the optimal solution for the IPl 
cannot be more than n. Hence the integrality gap is Q(n). □ 

High integrality gap necessitates an alternative approach. 

Approximate Integer Program We propose an alternative integer program 
IP 2 which is a form of Lagrangian relaxation of IPl. We will show that its 
maximal solutions are also solutions of IPl and any a approximation of IP2 is 
a 2a/(l — e) approximation of IPl. A maximal solution of IP2 is a solution in 
which z v — max{0, Y^eeSM V&~ c v\ for all v and deletion or addition of an edge 
does not improve the objective function value. 

IP2: maximize </> = 2 y e — (1 + e) z v , for some e > 

eeE vev 

subject to Ve < 

z v G {0,1,2,...} Vv G V 
y e G {0,1} Ve G E 

Note that any subset of edges E' is a feasible solution of IP2 if we choose z v = 
max{0, YleeSfv) V e ~ Cl >) ^ or au v - Besides these values of z will give maximum 
value of the objective function. Hence z values are not required to be specified 
in the solutions of IP2. 



IPl: maximize 
subject to 



Lemma 4. Every maximal solution of the integer program IP2 is also a feasible 
solution of PDBEP. 

Proof. Consider any maximal solution E' of IP2. In a maximal solution z v = 
max{0, X)ee6M V e ~ Cy } ^ or a ^ v - Assume that it is not a feasible solution of 
PDBEP. Then there must exist an edge e = (u, v) £ E 1 such that z u > 1 and 
z v > 1. Define an alternative solution E" = E' \{e} and decrement z u and z v 
by 1 each. Observe that the objective function of the new solution increases by 
2e. This contradicts that E' is a maximal solution. □ 

Lemma 5. Any a approximate solution of IP2, which is also maximal, is a 
2a/ (1 — e) approximation of PDBEP problem. 

Proof. Let E' be an a-approximation maximal solution of IP2 with m% — z' v 
and mi = J2 e y'e ~ TO 2 = \E'\ — m-2. Let E" be an optimal solution of PDBEP. 
Then y" = 1 for e £ E" and z" = max{0, J2e£8(v) v'e~ c A is a solution of IP2, i.e., 
E" is also a solution of IP2. Define n 2 = X]„ 2 " an d n i = S e 2/" ~ n 2 = ri2- 
We have = 2mi + (1 - e)m 2 and </>(£") = 2rii + (1 - e)n 2 . Let OPT 

denote the optimal value of the IP2 objective function. Then <fi(E") < OPT and 
OPT/a < cf>(E'). So 2ni + (1 - e)n 2 < a(2mi + (1 - e)m 2 ). So (1 - e)\E"\ = 
(l-e)(ni+n 2 ) < 2ni + (l-e)n 2 < a(2mi+(l-e)m 2 ) < 2a(mi+TO 2 ) = 2a-\E'\. 
□ 

2.4 Algorithm for IP2 

We propose Algorithm [3] which approximates the IP2 problem within a constant 
factor of approximation. LP2 is the linear program relaxation of IP2. Here we 
assume that an additional constraint is imposed, namely, {z v = 0\v £ C} where 
we require a solution in which every v € C must necessarily satisfy the degree 
constraint. The input to the problem is (H = (V, E),C). Algorithm starts with 
E' = and builds it up one edge at a time by iterative rounding. In each iteration 
we discard at least one edge from further consideration. Hence it requires at 
most \E\ iterations (actually it requires at most |^| + 1 iterations, see the remark 
below.) In the interest of ease in analysis Algorithm[3]is presented in the recursive 
format. 



LP 2: maximize 




z v , for some e > 



subject to 





z v >0WveV 
y e > Ve £ E 
-y e > -1 Ve £ E 



Data: A connected graph G = (V, E) and a function c : V — > N 
Result: A solution of PDBEP problem, 
for v £ V do 

I JV - = C V ] 

end 

C:=0; 

E' := SolveIP2{G,C,f); 
return E'; 

Function: SolveIP2(# = (V H ,E H ),C,f) 
if E H ~ then 

| return 0; 
end 

Delete all isolated vertices from Vh', 
(y,z) = LPSolver(#,C); 

/* solve LP2 with degree-bounds f(x) for all x € Vh */ 

if 3e £ Eh with y E = then 
Hi := (^H,i5 H \{ e }); 
Ci := C; 

£' := SolveIP2{H 1 ,C 1 J); 
else 

From Lemma[6]there exists an edge e := (u,v) with y e > 1/2; 
From Lemma [7|w.l.g. we assume (/„ > 0, z v = 0); 

fv := /tJ — lj 

Ci -CUW; 

/„ := max{/„ - 1,0}; 

ffi := \{e}); 

E' := SolveIP2(Hi, d,f)U {e}; 

/* By including e in E we effectively rounded up y e to 1. Hence 
in case f u = then implicitly z u is also raised to ensure 
that the condition X^ e 'g<5(u) fe' ^ fu + z u continues to hold. We 
do not explicitly increase z u value since it is not output as 
a part of the solution. */ 

end 

return E'\ 

Algorithm 3: Iterative Rounding based Algorithm in Recursive Format 



In the following analysis we will focus on two problems: (H, C) of some i-th 
nested recursive call and (ifi, C\) of the next call. For simplicity we will refer to 
them as the problems of graphs H and H\ respectively. 

Lemma 6. In a corner solution of LP2 on a non-empty graph there is at least 
one edge e with y e = or y e > 1/2. 

Proof. Assume the contrary that in an extreme point solution of LP2 all y e are 
in the open interval (0,1/2). Let us partition the vertices as follows. Let n\ 
vertices have f v > and z v > 0, n 2 vertices have f v > and z v — and n3 
vertices have f v = Q and z v > 0. Note that the case of /„ = and z v = cannot 
arise because y e > for all e. In each case let n\ vertices have the condition 
Se65(u) He < fv + %v tight (an equality) and n" vertices have the condition a 
strict inequality. Let the number of edges be m. 

The total number of variables is ni+n 2 +n 3 +m. In n^+n^ cases SeeSM ^ e = 
f v + z v where f v >l and each y e < 0.5 so there must be at least 3 edges incident 
on such vertices. Since the graph has no isolated vertices, every remaining vertex 
has at least one incident edge. Hence m > (3^^ + 3n' 2 + n'{ + n 2 + n 3 )/2. So 
number of variables is at least n[ + n' 2 + (1.5)(rii + n 2 + n^). 

Now we find the number of tight conditions. None of the y e touch their 
bounds. The number of z v which are equal to zero is n 2 , and the number of 
instances when J2ee8(v) = fv + z v is n\ + n' 2 + n' 3 . Hence the total number 
of conditions which are tight is n 2 + n\ + n' 2 + ji' 3 . Since the solution is an 
extreme point, the number of tight conditions must not be less than the number 
of variables. So n 2 + n[ +n' 2 +n' 3 > n[ +n' 2 + (1.5)(ni + n i + n z)- This implies 
that ni = n 2 = n 3 =0, which is absurd since the input graph is not empty. □ 

Remark: The program LP2 has \E\ + \ V\ variables and 2\E\ + 2\V\ constraints. 
Hence in the first iteration the optimal solution must have at least \E\ — \V\ 
tight edge-constraints (i.e., y e = or y e = 1.) All these can be processed simul- 
taneously so in the second round at most \V\ edges will remain in the residual 
graph. Thus the total number of iterations cannot exceed |V| + 1. 

Lemma 7. If y e > 1/2 in the solution of LP2 where e — (u,v), then (f u > 
0,z u = 0) or (f v >0,z v = 0). 

Proof. Assume that z v > and z u > in the solution. Let minimum of z v , z u , 
and y e be j3. Subtracting (3 from these variables results in a feasible solution 
with objective function value greater than the optimum by 2/3 • e. This is absurd. 
Hence z u and z v both cannot be positive. 

Next assume that f u — and z u — 0. Then y e must be zero, contradicting 
the fact that y e > 1/2. Similarly /„ = and z v = is also not possible. 

Therefore either (/„ > 0, z u = 0) or (/„ > 0, z v = 0). □ 

Lemma 8. The Algorithm^ returns a feasible solution of PDBEP. 

Proof. The claim is trivially true when the graph is empty. We will use induction. 



In the case of y e — 0, the solution of H\ is also the solution of H. From 
induction hypothesis it is feasible for H\ hence it is also feasible for H. 

Consider the second case, i.e., y e > 1/2. Let e = (u,v). From Lemma 
= a > and z v — 0. Since z\ v = and f\ v = a — 1, in the solution of Hi at 
most a — 1 edges can be incident on v. So there are at most a edges incident on 
v and f v = a in the solution of H . Thus e is valid in the solution of H. Other 
edges are valid due to induction hypothesis. □ 

Now we analyze the performance of the algorithm. 

Lemma 9. Algorithm^ gives a 1.5/(1 — e) approximation of IP2. 

Proof. Let c denote 1.5/(1 — e). We will denote the optimal LP2 solutions of H 
and Hi by F and F\ respectively. Similarly I and Ii will denote the solutions 
computed by the algorithm for H and Hi respectively, /i* and zi* denote the 
parameters associated with Hi, We will assume that z x = max{0, Y2ee5(x) V e ~ 
f x } for integral solutions to compute their 0- values. Again we will prove the 
claim by induction. The base case is trivially true. From induction hypothesis 
<f>(Fi) I '(f)(h) < c and our goal is to show the same bound holds for (f>(F) / (f>(I) . 

In the event of y e = in F, (f>(F) = </>(F±) and (f)(1) — (f>(h). Hence 
(f>(F)/(f>(I) = (f>(Fi)/(f)(I 1 ). 

In case y e = a > 1/2 we will consider two cases: (i) f u > and (ii) /„ = in 
F. In the first case / differs from Ii in three aspects: y e = 1 in /, f v = fi v + 1 
and f u — fi u + 1. So z v and z u remain unchanged, i.e., z v = z\ v and z u = zi u . 
Thus (f>(I) — (f)(Ii) + 2. In the second case also y e increases by 1 and z v remains 
unchanged but z u increases by 1 because in this case f u = fi u — 0. Hence 

(f)(1) =(f)(Il) + 1-6. 

From induction hypothesis <f>(Fi) < c<fr(Ii). Hence for any solution F{ of Hi, 
we have (f>(F{) < c<p(Ii). In the remaining part of the proof we will construct a 
solution of LP2 for Hi from F, the optimal solution of LP2 for H. 

Again we will consider the two cases separately. First the case of /„ > 0. Set 
y e = 0. If X) e 'e5(i))\{e} y^' — 1 — a then subtract the values of y e > fore' £ 5(v)\{e} 
in arbitrary manner so that the sum J2 e 'G5(v)\{e} V e ' decreases by 1 — a. If 
Se'e5(ii)\{e} 2/c < 1 — en, then set y e i to zero for all edges incident on v. Repeat 
this step for edges incident on u. Retain values of all other variables as in F 
(in particular, the values of z u and z v ). Observe that these values constitute a 
solution of LP2 for Hi . Call this solutionis'. Then <f>(F[) > <j>(F) -2(1 + 1 -a) > 
(f>(F) - 3. We have (f>(F) < <f>(F[) + 3 < c(0(I x )) + 3 < c(4>(I) - 2) + 3 < c(f>(I). 

In the second case f u = 0. Once again repeat the step described for edges 
incident on v and set y e to zero. In this case z u > a so subtract a from it. It is 
easy to see that again the resulting variable values form an LP2 solution of Hi, 
call itF[. So <j>(F[) = (f>(F) — (2 - (1 + e)a). So (p(F) = (j>(F[) + (2 - (1 + e)a) < 
c(f)(Ii) + (2 — (1 + e)a). Plugging (f>(I) — 1 + e for (f)(Ii) and simplifying the 
expression gives 4>(F) < c(f>(I). This completes the proof. □ 

Combining lemmas [5] and |H] we have the following result. 



Theorem 3. Algorithm^ approximates PDBEP with approximation factor 3/ (1 — 



3 Approximation Algorithm for the weighted case 

Let H(v) denote the heaviest c v edges incident on vertex v, called heavy set 
of vertex v. Then from a generalization of Lemma Q] the optimum solution of 
PDBEP in weighted-edge case is bounded by YlveV SeeffM u; ( e ) wnere w ( e ) 
denotes the weight of edge e. We will describe a method to construct upto 

1 + 2 log \V\ solutions, which cover L) v& vH(v). Then the heaviest solution gives 
a 2 + 2 log \ V\ approximation of the problem. 

3.1 The Algorithm 

Input: A graph (V, E) with non-negative edge-weight function w(). Let |V| = n. 
Step 1: E x = E \ {e = (u, v) e E\e £ H(u) and e £ H(v)}. 
Step 2: T = {e = (u,v) e E\e G H{u) and e € H(v)}. 

Step 3: E2 = E%\ T. Clearly each edge of E2 is in the heavy set of one of its 
end- vertices. Suppose e = (it, v) S E with e ^ H(u) and e G H(v). Then we will 
think of e as directed from u to v. 

Step 4: Arbitrarily label the vertices from to n — 1. Define sets of edges 
Ao, . . . , Ak~i and Bq, . . . , B^-i, where fc = logn, as follows. A r consists of edges 
(u,v) directed from u to v, such that the r — 1 least significant bits of binary 
expansion of the labels of u and v are same and r-th bit of u is zero and the 
same bit of v is one. B r is defined similarly except the r-th bit of u is one and 
that of v is zero. 

Step 5: Output that set among the 21ogn + 1 sets, T, Aq, . . . , B , . . . , Bk, 
which has maximum cumulative edge weight. 

Theorem 4. The algorithm gives a feasible solution with approximation factor 

2 + 2 logn. 

Proof. Set T constitutes a feasible solution since both ends of each edge in it 
satisfy the degree constraint. In A r all arrows are pointed from u with r-th 
bit zero to v with r-th bit one. Hence it is a bipartite graph where all arrows 
have heads in one set and the tails in the other. All vertices on the head side 
satisfy the degree conditions because all their incident edges are in their heavy 
sets. Therefore A r are feasible solutions. Similarly all B r are also feasible. It 
is easy to see that every edge of E2 belongs to A r or B r for some r. Hence 
TU {U r A r ) U (ll r B r ) = Ei. Observe that l) v H(v) = E x . Only T-edges have both 
ends in heavy sets. Using the fact that OPT < J2v w (H{v)), we deduce that 
OPT < 2w{T) + J2r( w ( A r) + w(B r )). So the weight of the set output in step 5 
is at least OPT/ (2 + 2 logn). □ 



4 Exact Algorithm 



In this section we give a polynomial time exact algorithm for the unweighted 
PDBEP problem for the special case when the input graph is a tree. We will 
denote the degree of a vertex v in the input graph by d(v) and its degree in a 
solution under consideration by d'(v). 

Let T be a rooted tree with root R. For any vertex v we denote the subtree 
rooted at v by T(v). Consider all feasible solutions of PDBEP of graph T(v) 
in which degree of v is at most c v — 1, call them iJ-solutions. Let h(v) be the 
number of edges in the largest such solution. Similarly let g(v) be the optimal 
G-solution in which the degree of v is restricted to be equal to c v . Lastly b(v) will 
denote the optimal i?-solution which are solutions of T(v) under the restriction 
that degree of v be at least c v and every neighbor of v in the solution satisfies 
the degree condition. It may be observed that one class of solutions of T(v) 
arc included in G-solutions as well as in _B-solutions. These are the solutions in 
which d'(v) = c v and every neighbor u of v in the solution has d'(u) < c u . If in 
any of these cases there are no feasible solutions, then the corresponding optimal 
value is assumed to be zero. Hence the optimum solution of PDBEP for T is the 
maximum of h(R), g(R), and b(R) and all three values are zero for a leaf nodes. 

Let Ch{v) denote the set of child-nodes of v in T(v). We partition Ch{v) 
into H(v) = {u e Ch{v)\h{u) > max{ 5 (u), b(u)}}, G(v) = {u e Ch(v)\g(u) > 
max{h(u),b{u)}}, B(v) = Ch(V) \ (G{v) U H(v)). While constructing a G- 
solution of T(v) from the solutions of the children of v we can include the edge 
(v, u) for any vertex u in H{v) U B(v) along with the optimal solution of T(u) 
without disturbing the degree conditions of the edges in this solution. But we 
can add edge (v,u) to the solution, for any u G G(v), only by selecting a B- 
solution or an iJ-solution of T{u) because if we use a G-solution for T(u), then 
vertex u which was earlier satisfying the degree condition, will now have degree 
Cu + 1- Similarly while constructing a i3-solution of T(v) we can connect v to 
any number of H(v) vertices and use their optimal //-solutions. But in order to 
connect v with u € B(v) U G(v) we must use the optimal if-solution of T(u). 

If k = c v - \H{v)\ - \B{v)\ > 0, then we define S'(v) to be the set of 
k members of G(v) with smallest values of g(u) — max{/i(u), b(u)}. Otherwise 
S'(v) = 0. Similarly if k = c v - \H(v)\ > 0, then we define S"(v) to be the set 
of k members of G(v) U B(v) with smallest key values where key is g{u) — h{u) 
for u e G(v) and b(u) - h(u) for u E B{v). Otherwise S"(v) = 0. Now we have 
following lemma which leads to a simple dynamic program for PDBEP. 

Lemma 10. For any internal vertex v ofT, 

(i) h(v) = J2ueB(v) b (u)+J2ueH(v) h ( u ) + J2ueG(v) fl(w)+min{c„-l, \H(v)\ + 
\B(v)\}, 

If d(v) — c v and v ^ R, then set b(v) = g(v) = otherwise 
(H) 9{v) = EueB(v) b (u) + EueH(v) h (u) + E„ 6 G(»)\S'(») + 
J2ues'(v) max{^(u), b(u)} + c v , 



(Hi) b{v) = T,ueH( v )us"(v) h (u) + Y,ueB(v)\s"(v) h ( u ) +T, u eG(v)\s»(v) 9(u) + 
max{c„, \H(v)\}, 

Observe that if h[u) is equal to b(u) or g(u), then u is categorized as an H (v) 
vertex and if b(u) = g(u) > h(u), then u is assigned to B(v) set. Hence the last 
term is maximum in each of the cases in the lemma. 

The algorithm initializes h(v),b(v), and g(v) to zero for the leaf nodes and 
computes these values for the internal nodes bottom up. Finally it outputs the 
maximum of the three values of the root R. Computations for any internal vertex 
takes 0(|C/i| log \ Ch\) time where Ch is the set of children of that vertex. Besides 
the ordering the vertices so that child occurs before the parent (topological sort) 
takes 0(n) time. Hence the time complexity is O(nlogn). 

5 Future work 

It remains an open question if there exists a constant factor approximation algo- 
rithm for the PDBEP problem when the input graph is weighted. The objective 
function of LP1 can be easily modified to handle the weighted case but due 
to the large integrality gap this approach remains useless. However, there are 
cutting-plane methods like Chvatal-Gomory cuts [J] that have been known to 
improve the integrality gaps for some problems. It would be worthwhile to see if 
these methods can help reduce the integrality gap of LP1. 

There seems no reason to assume that the number of resources necessary for 
the accomplishment of a job in the resource allocation problem cannot exceed 
two. Hence a natural generalization of the PDBEP problem to hypergraphs. 

Any NP-hard problem is not considered resolved unless and until it has an 
algorithm with approximation factor that matches the lower bound for that 
problem. As far as we know there is no known inapproximability result for the 
PDBEP problem. So that presents another avenue for further research. 
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